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Recent observation of zero bias conductance peaks in semiconductor wire/superconductor heterostructures 
has generated great interest, and there is a hot debate on whether the observation is associated with Majorana 
fermions (MPs). Here we study the local and crossed Andreev reflections in a junction of two normal leads and 
a sandwiched superconductor-semiconductor wire with two spatially separated but strongly coupled MF end 
states. The conductance and Fano factors of such a device are sharply different from the topologically trivial 
case even in the presence of disorder, and can hence be used to identify MFs unambiguously. 
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Introduction — Recent proposals suggest that Majorana 
fermions (MFs) can appear as zero energy end states in su- 
perconductors constructed by inducing superconductivity on 
semiconductor wires with Rashba spin-orbit coupling 

Remarkably, several experimental groups [8-I0I re- 
cently reported the observation of zero bias conductance 
peaks (ZBCPs) in Andreev reflection experiments from a 
normal lead to the end of the aforementioned semiconduc- 
tor/superconductor heterostructure. These ZBCPs were possi- 
bly due to the MF induced Andreev reflection fTT','T2']. How- 
ever, the oiigin of these ZBCPs remains a subject of debate 

iBlil. 

In this work, instead of studying MF induced local Andreev 
reflection as done in previous works fo'-lo'], we explore the 
non-local properties of MFs. Since each MF has only half the 
degrees of freedom of a Dirac fermion and two spatially sepa- 
rated MFs can form a non-local Dirac fermion, such non-local 
Dirac fermions can induce non-trivial cuiTent-current correla- 
tions between two spatially separated normal leads which are 
coupled to the two MFs. To be specific, an experimental setup 
depicted in Fig. la is studied, in which two normal leads are 
attached to the two ends of an s-wave superconductor with 
Rashba spin-orbit coupling in the presence of a magnetic field 
parallel to the wire. 
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FIG. 1 : (a) A schematic setup of experiment, two normal leads are 
coupled to the two ends of a superconductor which supports MFs. 
(b) A local Andreev reflection process. An electron from one lead is 
reflected as a hole in the same lead, (c) A crossed Andreev reflection 
process. An electron from one lead is reflected as a hole in another 
lead and a Cooper pair is injected into the superconductor. 



When the wire is in the topological regime with MF end 
states and if the MF end states from the two ends of the wire 
are coupled, local Andreev reflection processes (as depicted 
in Fig. lb) can be suppressed and the MF end states induce 
crossed Andreev reflections, in which an electron from one 
lead is reflected as a hole in a different lead (as depicted in 
Fig.lc). Therefore, each normal lead tunnels one electron into 
the superconductor in each tunnelling event. As a result, the 
Fano factor, which is the ratio of the shot noise of the cur- 
rent and the average tunnelling cuiTent, is e in the topologi- 
cal regime. On the other hand, the Fano factor which mea- 
sures the charge transfer of each tunnelling event, is 2e in 
the topologically trivial regime when local Andreev reflection 
processes dominate. Therefore, the measurement of the shot 
noise can be used to distinguish between the topological and 
non-topological regimes of the wire unambiguously. 

It is important to note that an experimental setup similar to 
Fig. 1 a has been fabricated recently 

mm 

. The shot noise of 

leads of the set up is measured f20'l . However, the experiment 
focused only on the trivial regime with zero external magnetic 
field and the strong magnetic field regime where the proximity 
gap is suppressed. 

Model and Formalism — To model the quasi-one dimen- 
sional s-wave superconductor with Rashba spin-orbit coupling 
as shown in Fig. la, we use the following tight-binding model 
ifTi flsIl with Nx and Ny sites in the x and y directions respec- 
tively: 

HqlD = ^ -<(V^R+d,a^fl," + h.C.) - llTpl^^^^pR.^ 

+ I^R,d,a./3 ^^^R'^R+'i,J ■ ^ d)„^^R,^ + h.C. 

Here, R denotes the lattice sites, d denotes the two unit vec- 
tors dx and dy which connects the nearest neighbor sites in 
the X and y directions respectively, a, f3 are the spin indexes. 
t is the hopping amplitude, /i is the chemical potential, Uu is 
the Rashba coupling strength, Vx is the Zeeman energy caused 
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by a magnetic field along the wire in the x-direction. A is the 
superconducting pairing amplitude and Vimp(R) is the on-site 
random impurity which is randomly distributed in the range 
[— w/2, w/2] at each site. In this work, we set = 2A such 
that the superconducting wire can support MF end states by 
tuning the chemical potential. 

The parameters in the tight-binding model are chosen to 
match the corresponding values in a recent experiment done in 
Ref-lUl. Here, A = 250^eF, t = 25 A and Ur = 2.5 A. The 



dimensions of the wire are N^a « l/im and NyU w lOOnm 
which is about twice the superconducting coherence length 
of the system. The length of the wire chosen is about half 
the length of the experimental value in Ref.[8]. Due to the 
short length of the wire, as shown in Fig. 2a, the energy versus 
chemical potential plot exhibits oscillatory behaviour in the 
topologically non-trivial regime |19, 26] instead of staying at 
zero energy as the two MF end states can couple to each other 

To study the current-current correlation mediated by the 
MF end states, two semi-infinite normal leads are attached to 
the two ends of the TS as shown in Fig. la. The two normal 
leads are described by Eq{T|by setting A to zero. The tun- 
nelling barriers are simulated by the reduced hopping ampli- 
tudes iic = tuc — 0.3t between the leads and the supercon- 
ductor where t^c (tnc) denotes the hopping amplitude from 
the left (right) lead to the TS. 

We use the recursive Green's function method to calculate 
the scattering matrix of the model |21] where the scattering 
matrix is related to the Green's functions by 
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Here, S'^j'^ is an element of the scattering matrix which de- 
notes the scattering amplitude of a /3 particle from lead j to an 
a particle in lead i. i,j — 1 or 2. 1 and 2 denote the left and 
the right lead respectively, a, /3, G {e, h} denotes the electron 
(e) or hole (h) channels. G"" is the retarded Green's function 
of the superconducting wire. Tf — i[(Sf )'' — (Ef )°], where 
(j]ay{a) jjjg p^ficjg retarded (advanced) self-energy of 
lead i. 

With the scattering matrix, the average current /, the dif- 
ferential shot noise Pij and shot noise can be calculated 



j[r,"]i/2=,G'-*[rf]i/2. 
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as: 1221, 121 

h = f /;''ETr[/-sgn(a)^-(£;)t^-(£;)]d£;, 



p, 



E sgn(a)Tr[5:f {EY St,' iE)Sf,' Sf,' {E)l 
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where sgn(a) = 1 if a = e and sgn(a) = — 1 if a = h . 
In this work, we set the chemical potential of the two nor- 
mal leads to be the same and the voltage bias between the 
leads and the superconductor to be V. Physically, Cy = 
J-^ {0)SIj {t)dt measures the current fluctuation of leads 
i and j, where 5Ii = Ii{t) — It denotes the deviation of the 
current at time t with respect to the average current Ii. At 
low temperatures with T <C eV/ks, the current fluctuation is 
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FIG. 2: (a)The energy spectrum of a quasi-lD wire with Nj:a = 
l/im. The topological region is indicated above, (b) Contour plot of 
differential conductance Gi of the left lead as a function of chemical 
potential and electron incident energy E. The conductance shows an 
oscillatory behavior in the topological regime, (c) Gi versus incident 
energy at a fixed chemical potential denoted by the red dashed line 

in (b). The height of the peak at E ^ Em is 



dominated by the shot noise |24t] and thus we simply denote 
Cij as the shot noise and Pij as the differential shot noise. 

When the superconducting wire is in the topological regime 
with only one transverse subband occupied and two MFs ap- 
pear as end states, we expect the transport properties of the 
N/TS/N junction at eV <C A to be qualitatively described by 
the effective Hamiltonian H^f / = Hl + Hm + Ht, where 

Hm = iEMlLlR 

Ht = -i[tLlL{^l{0) + V'l(O)) + ifl7fl(^l?(0) + V'fl(O))]. 

(4) 

In EqlH ipL/R denotes a fermion operator of the left (right) 
normal lead, v f is the corresponding Fermi velocity. Em is 
the coupling strength between the two MF end states and 
•yji. The coupling strengths between the normal leads and the 
TS are denoted by and Ir respectively. This model was 
first introduced in Refs. ll24 12511 and the scattering matrix of 
the Hamiltonian can be found using the equation of motion 
approach [11]. With the scattering matrix, the current and shot 
noise can be calculated using EqlS] In the following, we show 
that the numerical results of the tight-binding model in the first 
topological regime, where only one transverse subband of the 
wire is occupied, and the results from the effective Hamilto- 
nian model agree well with each other However, the effective 
model Hef / does not apply to the wires with more than one 
occupied transverse subband. The multi-channel regime and 
the effect of disorder are studied based on the lattice Green's 
function approach. 

Current — In this section, we focus on the tunnelling 
current near the topological regime where only one trans- 
verse sub-band is occupied. However, due to the short 
wire geometry, the MFs at the two ends couple to each 
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other quite strongly. Due to the oscillatory nature of the 
MF wavefunctions, the coupling strength of the MFs oscil- 
lates and the resulting coupling energy oscillates as Em « 
h kp ^-2N^a/io cosikpNx), whcrc kp is the Fermi momen- 
turn which is a function of chemical potential and magnetic 
field lE6l] . The energy spectrum of the TS is shown in Fig. 2a. 
The topological regime marked out in Fig. 2a is determined in 
the long wire limit where zero energy modes appear As the 
chemical potential increases and a second transverse subband 
is occupied, the TS becomes topologically trivial. 

To study the MF end states, we calculate the differential 
conductance dli / dV using Eql3] The contour plot of the dif- 
ferential conductance Gi — dli/dV as a function of elec- 
tron incident energy E = eV and the chemical potential 
is shown in Fig. 2b. As expected, the MFs manifest them- 
selves by inducing conductance peaks. However, when the 
two MFs couple to each other and Em ^ ^1,^2 ' the height 
of the differential conductance peak is reduced to Gi {Em) = 
dIi/dV\ev=EM ~ 2e 1^ from the usual 
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lead geometry. Here, ti = y-^^L ih = y-^tn) is the ef- 
fective coupling strength between the left (right) normal lead 
and the TS. 

Another interesting point for tunnelling into a TS with two 

strongly coupled MFs is that at low incident energy E -C Em, 

2 ^2^2 

the differential conductance Gi (E) w ^ -kr- depends on the 
product of ti and ^2- This means that an electron from a nor- 
mal lead cannot tunnel into the TS unless a second normal 
lead is coupled to the TS. This is a manifestation of the fact 
that local Andreev reflection processes are suppressed and the 
current is purely caused by crossed Andreev reflection pro- 
cesses in this regime. In other words, the strongly coupled MF 
end states drastically change the local and non-local transport 
properties of the N/TS/N junction. As it is shown below, the 
measurements of the differential shot noise and shot noise can 
be used to probe the topological regime unambiguously. 

Differential Shot Noise — The contour plots of the dif- 
ferential shot noise Pii{E) and Pi2{E) defined in Eql3]are 
shown in Fig. 3a and Fig. 3b, respectively. In Fig. 3, Pu 
and P12 at energy E are normalized by the differential con- 
ductance Gi{E) and the average differential conductance 
G{E) = i(Gi + G2) respectively. The ratio F{E) = 
Pii{E) I Gi {E) is the Fano factor for electrons at incident en- 
ergy E. At weak tunnelling limit, the Fano factor gives the 
value of the tunnelling charge at each tunnelling event. From 
Fig. 3a and Fig.3e, it is clear that the Fano factor at i? = is 
the electron charge e in the topological regime. Similarly, the 
Fano factor for the right lead is P22 {E = 0) /G2(£; = 0) = e. 
This indicates that for each tunnelling event, each normal lead 
contributes one electron in the tunnelling process. Moreover, 
it is evident from Fig. 3b and Fig.3e that Pi2{E = 0)/G = e, 
where P12 is the current-current correlation of two differ- 
ent leads. As pointed out in Ref.[24], the cross correlator 
P12 is bound by the relation 2|Pi2| < Pn + P22 for any 
stochastic processes. At E — Q, we have 2|Pi2(i?)| = 
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FIG. 3: (a) Contour plot of Fano factor Pn/eGi for electrons 
with incident energy E at chemical potential /i. (b) Contour plot 
of the crossed Fano factor P12 /eG. (c) Contour plot of Fano factor 
Cii/e/i as a function of voltage bias eV and chemical potential, (d) 
Contour plot of crossed Fano factor C12 /el. (e) The Fano factors 
Pii/eGi and Pi2/eG as a function of incident energy i5 at a fixed 
chemical potential denoted by the dashed lines in (a) and (b). (f) The 
Fano factors Cn/e/i and Ci2/eJ as a function of voltage bias at 
fixed chemical potential denoted by the dashed lines in (c) and (d). 

Pll{E) + P22{E). 

This indicates that the two leads are perfectly correlated 
with each other such that a Cooper pair is injected into the su- 
perconductor at each tunnelling event (as depicted in Fig.lc). 
This is in sharp contrast to the topologically trivial regime as 
shown in Fig. 3a in which local Andreev reflection processes 
(as depicted in Fig. lb) dominate and the Fano factor for each 
lead is 2e instead. The tunnelling currents of the two leads are 
only weakly correlated in the absence of MFs. 

Away from P = 0, we see that the Fano factor in the topo- 
logical regime is strongly energy dependent. It is also shown 
in Fig. 3a and Fig.3e that the Fano factor Pii(P)/Gi(P) 

(2 

reaches a minimum value of .■i?..^ e at P = Em when the 
tunnelling conductance is maximum. For P ^ Em, local 
Andreev reflection processes dominate and the Fano factor 
approaches 2e. On the contrary, the shot noise of 2e in the 
trivial regime is insensitive to the energy of the incident elec- 
trons as long as the energy is much smaller than the first An- 
dreev bound state energy of the wire. Since the Fano factor 
behaves very differently in the topologically trivial and non- 
trivial regimes, it can be used to determine the topological 
regime unambiguously. As it is shown in Fig. 2a, Em is about 
0.1 A for a realistic wire and it can be even larger for a shorter 
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wire Therefore, the experimental measurements of the 
energy dependence of the Fano factor is highly feasible. 

It is important to note that the cross correlator P12 also ex- 
hibits strong energy dependence in the topological regime. As 
shown in Fig. 3b and Fig.3e, P^/eG « "a, a ~ 1 for 

E < Em- AtE^ Em, Pii/eG « -j^^- At £; > Em, 
local Andreev reflections dominate at each lead and the corre- 
lations between the two leads drop to zero eventually. On the 
other hand, in the absence of MFs outside of the topological 
regime, P12 drops to zero quickly, independent of the incident 
energy. Therefore, the measurement of P12 provides another 
way to probe the topologically non-trivial regime. 

Shot Noise — In this section, we study the shot noise dj , 
which is the integration of the differential shot noise over 
the incident energy of the electrons as defined in Eq|3] The 
contour plots of Cn and C12, normalized by Ii and / — 
+ I2) respectively, as a function of chemical potential 
and voltage bias are shown in Fig. 3c and Fig. 3d. The Fano 
factor Cii I Ii gives the charge leaving lead i at each tunnelling 
event. As expected, in the crossed Andreev reflection regime 
with E <C Em, Ca/Ii — e. In this case, the two leads are 
perfectly correlated as C12/I = e. 

Outside of the topological regime, Cu/Ii ~ 2e as local 
Andreev reflection processes dominate. Moreover, the cross 
correlation between the two leads C12 is significant only when 
the incident energy of electrons satisfies E <^ Em in the topo- 
logical regime as shown in Fig. 3d. 

Disorder effect — The observation of the ZBCPs in recent 
tunnelling experiments is an important step in the search for 
MFs in solid state systems fS MlOll . However, as argued in 
Ref.|13], fermionic end states formed by two MF end states 
in the topologically trivial regime can also induce ZBCPs in 
the presence of disorder Therefore, it is important to distin- 
guish the fermionic end states from the true MF end states. In 
this section, we show that the shot noise caused by a local- 
ized fermionic end state and the shot noise caused by a non- 
local fermionic state formed by two spatially separated MF 
end states are different. 

On-site random disorder with strength uj = 20A is added 
to the superconducting wire. The contour plots of the conduc- 
tance Gi and G2 are shown in Fig.4a and Fig.4b respectively. 
It is important to note that in the topologically trivial regime 
where two transverse subbands of the wire are occupied, a 
fermionic end state which has energy close to zero is induced 
by disorder at /i « — 84A. The ground state wavefunction 
at /X = — 84.6A is shown in Fig.4f and it is evident that the 
ground state is localized at the left end of the wire. As ex- 
pected, this zero energy fermionic end state induces a strong 
conductance peak for the left normal lead-superconductor 
junction. Moreover, the energy of the fermionic state oscil- 
lates as a function of chemical potential and it is difficult to 
distinguish this fermionic state from a true MF end state by 
measuring the conductance alone. 

However, since the fermionic state at « —84 A is a lo- 
calized state, the cross current-current correlation P12/G in- 
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FIG. 4: Random disorder with u = 20A is present for all the fig- 
ures, (a) and (b): The contour plots of the conductance Gi and G2, 
respectively, (c) and (d): The Fano factors Pn/e/i and Pvijel, re- 
spectively, (e) The ground state wavefunction magnitude I^'r^ in 
the topological regime with \i = —92. 7 A (indicated by dashed line 
A in (a)). The dimensions of the wire is A^^ = 50a and Ny = 5a. (f) 
The ground state wavefunction magnitude in the trivial regime with 
/i — — 84.6A (indicated by dashed line B in (a)). 

duced by this state is small as shown in Fig.4d. On the con- 
trary, P12/G is close to 1 in the topological regime. Moreover, 
the Fano facotor Pn/eGi at « is close to 1 only in the 
topological regime which indicates that crossed Andreev re- 
flection processes dominate. The ground state wavefunction 
in the topological reigme at /i = — 92.7A is shown in Fig.4e. 
It is evident that this fermionic end state, which can medi- 
ate crossed Andreev reflections, is a non-local fermionic state 
formed by two MF end states. 

Conclusion — We show that the measurements of current 
and shot noise in a N/TS/N junction can be used to determine 
the topological regimes of the superconducting wire with MF 
end states unambiguously. Particularly, the shot noise can be 
used to distinguish the fermionic end states from true MF end 
states even in the presence of disorder 
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